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Abstract 

^This  thesis  examined  the  electromagnetic  scattering  of  a  dihedral 
corner  reflector  for  both  monostatic  and  bistatic  source-receiver  geometries 
The  full  polarization  characteristics  of  the  electromagnetic  scattering  was 
maintained  while  the  dihedral  comer  reflector  was  modeled  on  and  off  the 
principal  (horizontal  reference)  plane.  The  dihedral  comer  reflector  was 
composed  of  two  perfectly  conducting  9  inch  square  flat  plates  joined  along 
a  common  edge  so  that  they  formed  an  interior  angle  of  90  degrees.  The 
frequencies  of  interest  included  2,  3,  and  10  GHz.  The  theoretical  modeling 
utilized  the  Uniform  Theory  of  Diffraction  4UTD>  and  the  Method  of 
Moments  (MOM).  The  analysis  was  performed  on  a  Cray  X-MP 
supercomputer  and  a  VAX  11/785  minicomputer.  The  results  from  the 
models  were  compared  with  each  other,  and  with  measurements  for 
valiu£i;on.  The  measurements  were  taken  on  the  principal  plane  for  same 
sense  and  cross  polarization  components.  The  MOM  model  compared  well 
with  measurement  for  the  monostatic  and  bistatic  same  sense  polarization 
components.  However,  differences  were  noted  for  the  cross  polarization 
component  comparison.  The  MOM  model  was  then  extended  to  off 
principal  plane  geometries,  which  was  not  meamred.  A  substantial 
amount  of  cross  polarization  scattering  was  predicted  for  off  principal 
plane  geometries. 
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FULL  POLARIMETRIC  NON-SPECULAR  SCATTERING 


OF  FLAT  PLATE  STRUCTURES 


Since  the  advent  of  radar,  engineers  and  scientists  in 


electromagnetics  have  attempted  to  predict  the  radar  echo  area  of  various 
targets.  The  Radar  Cross  Section  (RCS)  of  a  target  relates  this  echo  area  to 


the  radar.  In  order  to  predict  the  RCS  of  a  target,  one  needs  a  thorough 
understanding  of  the  scattering  principles  of  electromagnetic  energy.  The 
basis  of  all  electromagnetic  theory  begins  with  Maxwell’s  equations. 
However,  closed  form  solutions  are  available  only  for  simple  geometries. 
Attempts  to  extend  solutions  to  more  complex  scattering  problems  have 
brought  several  theories,  such  as  geometrical  optics,  physical  optics,  and 


the  theory  of  diffraction,  all  of  which  are  based  on  some  approximation. 


Present  radar  systems  operate  monostatically  (colocated  transmit 
and  receive  antennas)  and  use  same  sense  polarization  (same  receive  and 


transmit  polarization).  This  has  influenced  research  to  attempt  to  predict 


the  RCS  of  targets  for  monostatic  geometries,  and  same  sense  polarization. 
When  attempting  to  model  RCS,  a  target  is  decomposed  into  simple  shapes 
such  as  flat  plates,  cones,  cylinders,  etc.  The  scattering  from  each  of  the 
components  is  predicted,  with  the  complex,  vector  sum  yielding  the  total 
RCS  of  the  target.  The  scattering  from  simple  shapes  is  understood  for 


monostatic  geometries,  for  same  sense  polarization,  and  for  aspect  angles 
near  the  specular  (point  of  mirror-like  reflection). 

Future  radars  will  operate  both  monostatically  and  bistatically,  and 
will  likely  be  polarization  diverse.  These  requirements  dictate  that  the 
modeling  techniques  for  simple  shapes  be  extended.  These  extensions  will 
be  in  the  area  of  diffraction  terms  normally  insignificant  in  geometries 
and  polarizations  of  current  radar  systems. 

Background 

To  better  prepare  the  reader  to  understand  the  scattering  from 
corner  reflectors,  a  brief  historical  survey  is  in  order.  Comer  reflectors 
are  commonly  used  for  monostatic  radar  enhancement  (Kennaugh, 1956b). 
In  this  application,  echo  areas  much  larger  than  the  maximum  cross 
sectional  area  of  the  reflector  may  be  obtained. 

A  dihedral  (double  bounce)  corner  reflector  is  composed  of  two  flat 
plates  connected  to  form  a  comer,  thus  the  name  comer  reflector.  Double 
bounce  reflectors  do  not  reverse  the  rotational  sense  of  a  circularly 
polarized  wave,  and  are  therefore  desirable  for  use  with  circularly 
polarized  radars  (Hines, 1956a).  They  respond  to  circular  polarization  as 
well  as  to  linear  polarization  (Hines, 1956b). 

A  trihedral  (triple  bounce)  corner  reflector  is  composed  of  three  flat 
plates,  joined  to  form  a  comer.  These  operate  on  the  optical  principle  that 
an  arrangement  of  three  mutually  perpendicular  reflecting  planes  will 
return  a  triply-reflected  ray  in  the  direction  of  an  incident  ray 
(Baeumler,1958).  For  circularly  polarized  radars,  the  sense  of  the 
circularly  polarized  wave  is  reversed  upon  reflection  from  each  conducting 
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wall,  so  that  after  3  reflections,  the  returning  wave  is  of  opposite  sense  to 
the  incident  wave.  For  this  reason,  it  is  not  received  by  the  radar,  and  this 
is  why  trihedral  comer  reflectors  are  not  used  with  circularly  polarized 
radars  (Kennaugh, 1956a). 

There  are  several  components  in  analyzing  a  corner  reflector. 
There  are  specular  and  non-specular  radar  returns.  The  specular  return 
is  due  to  reflected  rays  which  obey  Snell's  Law  of  Reflection,  where  the 
angles  of  incidence  and  reflection  are  equal  with  respect  to  the  surface 
normal.  For  the  monostatic  radar,  the  reflected  ray  will  be  returned  in  the 
direction  of  the  incident  ray  (Baeumler,1958).  One  finds  that  the  solid 
angle  subtended  by  the  reflected  ray  is  inversely  proportional  to  the  square 
of  the  corner  leg  dimensions  expressed  in  wavelengths  (Kennaugh, 1956b). 
The  non-specular  return  is  due  to  diffracted  rays  which  result  from  edges 
and  surface  discontinuities.  The  diffracted  rays  will  scatter  in  many 
directions  and  usually  with  a  change  in  polarization. 

There  are  several  theories  attempting  to  explain  the  diffraction  of 
electromagnetic  energy.  The  Geometrical  Theory  of  Diffraction  (GTD) 
(Keller, 1962)  was  introduced  by  Keller  as  an  attempt  to  predict  diffraction. 
This  theory  had  its  share  of  shortcomings.  The  Uniform  Geometrical 
Theory  of  Diffraction  (UTD)  (Kouyoumjian,1974)  was  an  attempt  by 
Kouyoumjian  and  Pathak  to  overcome  these  shortcomings. 
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The  UTD  requires  keeping  track  of  all  the  diffracted  rays.  Because  of 
the  massive  bookkeeping  task  required  by  UTD,  most  solutions  include  only 
first  and  second  order  diffraction  terms,  which  correspond  to  singly  and 
doubly  diffracted  rays  respectively.  The  knowledge  of  solving  for  higher 
order  terms  is  available,  but  because  of  the  massive  bookkeeping  task 
required,  no  general  solution  includes  them.  (Marhefka,1984) 

Another  method  of  solving  RCS  problems  is  the  Method  of  Moments 
(MOM)  (Harrington, 1 966,1 968)  as  formulated  for  electromagnetic  field 
theory  problems  by  Harrington.  This  method  uses  numerical  techniques 
to  solve  Maxwell's  equations  in  the  integral  form,  but  the  electrical  size 
(size  in  wavelengths)  of  the  object  is  limited  by  computer  resources. 

Problem 

This  study  is  to  predict  the  non-specular,  electromagnetic  scattering 
of  flat  plate  structures  while  maintaining  polarization  diversity,  using 
modeling  techniques  which  handle  both  monostatic  and  bistatic 
geometries. 

Scope 

While  scattering  theories  are  valid  for  many  shapes,  this  study  will 
only  address  flat  plate  structures.  In  particular  the  dihedral  corner 
reflector,  which  is  a  scatterer  normally  found  on  targets  of  interest,  will  be 
analyzed  for  monostatic  and  bistatic  geometries,  aspect  angles  away  from 
the  specular  reflection  point,  and  cross  polarization.  Results  of  the 
analysis  will  be  critically  compared  with  measurements  of  the  same 
scatterer  to  validate  the  theoretical  development. 
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The  first  assumption  is  that  diffraction  terms  up  to  the  second  order 
are  all  that  are  needed  to  accurately  predict  the  scattering  from  the 
geometries  considered  in  this  study.  Second,  the  diffraction  coefficients 
used  in  this  study  are  assumed  to  be  accurate.  Finally,  the  MOM  is 
assumed  to  accurately  predict  cross  polarization  components. 

Standards 

The  experimental  data  obtained  by  The  Ohio  State  University 
ElectroScience  Laboratory  (OSU-ESL)  compact  range  will  be  used  as  the 
standard  to  which  the  theoretical  models  will  be  compared.  The  OSU-ESL 
compact  range  is  a  state  of  the  art  measurement  facility  which  can 
accurately  measure  electromagnetic  energy  from  2  to  18  gigahertz  (GHz). 
However,  this  research  will  require  bistatic  cross  polarization 
measurements  which  OSU-ESL  has  not  previously  measured. 

General  Approach 

The  problem  will  be  broken  into  steps  according  to  approach.  Step 
one  will  be  to  model  the  dihedral  (double  bounce)  corner  reflector  in  free 
space  using  the  UTD.  Step  two  will  be  to  model  the  dihedral  corner 
reflector  in  free  space  with  the  MOM.  Step  three  will  compare  the  UTD 
and  MOM  models  with  each  other.  Finally,  step  four  will  compare  the 
models  to  OSU-ESL  measurements 
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The  initial  attack  will  be  performed  with  the  use  of  a  two 
dimensional  UTD  corner  reflector  computer  code  developed  by  Aberegg 
and  Marhefka  at  OSU-ESL.  The  code  utilizes  geometrical  optics  along  with 
the  UTD.  The  code  includes  up  to  third  order  effects  which  include  single 
and  double  diffraction  terms. 

The  next  iteration  will  be  to  analyze  the  corner  reflector  with  the 
Electromagnetic  Surface  Patch  (ESP)  Code:  Version  II  -  Polygonal  Plates 
and  Wires  (Newman, 1985)  developed  by  Newman  at  OSU-ESL.  The  ESP 
code  allows  for  analysis  of  both  copolarization  and  cross  polarization.  The 
ESP  code  utilizes  the  MOM,  so  it  is  limited  in  use  by  the  electrical  size  of 
the  comer  reflector  in  conjunction  with  computer  resources. 

Materials  and  Equipment 

Hardware  needed  to  complete  this  study  includes:  a  VAX  11/785 
computer,  a  Cray  X-MP  computer,  and  a  graphics  plotter.  Software 
needed  includes:  the  ESP  code,  the  two  dimensional  UTD  comer  reflector 
code,  and  a  plotting  routine. 

Other  Support 

Support  from  Reed  and  Young  of  OSU-ESL  will  be  needed  in 
obtaining  experimental  data  to  be  used  for  comparison  with  theoretical 
results.  Support  for  the  use  of  the  computer  codes  will  be  provided  by 
Newman,  Marhefka,  and  Aberegg,  also  of  OSU-ESL.  Computer  support 
will  be  provided  by  the  Avionics  Laboratory. 
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A  survey  of  comer  reflector  RCS  prediction  is  presented  to  give  the 
reader  a  perspective  of  the  electromagnetic  scattering  for  this  geometry. 
In  Chapter  II,  general  concepts  of  electromagnetic  wave  theory,  diffraction 
theory,  and  the  moment  method  are  presented  in  relation  to  the  geometries 
of  interest.  Chapter  III  will  focus  on  the  analysis  of  the  scatterers  using 
the  UTD  and  MOM  codes.  Presentation  and  discussion  of  the 
experimental  analysis  and  comparison  to  measurements  will  be  covered  in 
Chapter  IV.  Finally  Chapter  V  will  present  conclusions  and 
recommendations. 

Summary 

The  purpose  of  this  study  is  to  predict  the  full  polarimetric,  non- 
specular  scattering  from  the  dihedral  corner  reflector  using  the  Uniform 
Geometrical  Theory  of  Diffraction  and  the  Method  of  Moments. 
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The  thrust  of  this  chapter  is  to  provide  the  reader  with  a  basic 
understanding  of:  1)  basic  theory  of  electromagnetic  waves;  2)  the 
Uniform  Geometrical  Theory  of  Diffraction  (UTD);  and  3)  the  Method  of 
Moments  (MOM).  The  discussion  will  be  limited  to  those  areas  applicable 
to  the  dihedral  comer  reflector. 

Throughout  this  writing,  ali  structures  are  assumed  to  be  perfect 
electric  conductors.  The  e  j®*  time  dependence  is  assumed  and  suppressed. 


Let  us  begin  this  discussion  with  a  definition  of  terms. 


=  Complex  Electric  Field  Intensity  Vector 
=  Complex  Scattering  Coefficient 

=  Horizontal  Unit  Vector 
=  Vertical  Unit  Vector 


where  the  bold  face  indicates  a  vector,  and  the  underbar  indicates  a 
complex  value,  and  the  "hat"  indicates  a  unit  vector. 
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An  arbitrarily  polarized  plane  wave  can  be  expressed  as  the  sum  of 
two  plane  waves  of  orthogonal  polarization: 

E  =  hEh  +  vEv  (2-1) 

where  the  subscripts  refer  to  the  respective  orthogonal  component,  which 
are  complex  in  general. 

The  radar  cross  section  (RCS)  of  an  object  is  an  equivalent  area 
defined  mathematically  by 

a  =  lim  4  7t  R2  ,ES|2/  |£i|2 

R^oo  (2-2) 


where 


a  =  Radar  Cross  Section 
R  =  Distance  Between  Radar  and  Object 
Es  =  Scattered  Electric  Field 
E1  =  Incident  Electric  Field 

Current  radar  systems  use  the  same  polarization  for  both  transmit 
and  receive  modes.  However,  the  RCS  is  polarization  sensitive  and  by 
decomposing  the  electromagnetic  field  into  the  orthogonal  components,  the 
RCS  can  be  decomposed  into  polarization  components.  The  scattered  field 
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can  be  related  to  the  incident  field  by  the  complex  scattering  coefficients  as 
follows: 


E^h  =  fihhE’h 

=  &hv  E’v 

ESvh  =  Syh  E’h 
Esw  =  SwE'v 


E^  =  ShhE’h  +  Shv  E’v 

Esv  =  £vh  E'h  "*■  Sw  E‘v 
In  matrix  notation  this  becomes 


- 

£hh  £hv 

E’h 

.Es,_ 

-  vh  &w  . 

LEV. 

The  polarization  scattering  matrix  is  now  defined  by 

' 

[  S  ]  =  ^hh  ^hv 
&vh  &w 


(2-4) 


(2-5) 


(2-6) 


and  the  RCS  is  related  to  the  scattering  matrix  by  (Knott, 1985) 

Oy  =  lim  4xR2  |fiij|2  (2-7) 

R-»oo 

where  the  double  subscript  refers  to  the  polarization  of  the  incident  and 
scattered  field,  i.e.,  i  =  h,v  and  j  =  h,v.  The  cross  polarized  terms  of  the 
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scattering  matrix  are  s^y  and  Sy},,  while  s^h  and  sw  are  the  copolarized 

terms.  The  target’s  horizontal  scattered  field  due  to  a  horizontal  incident 
field  is  defined  by  s^h.  and  the  target's  vertical  scattered  field  due  to  a 
vertical  incident  field  is  defined  by  svv.  The  target's  depolarization 
characteristics  are  defined  by  s^y  and  svh. 

The  scattering  coefficients  can  be  expressed  in  terms  of  the  incident 
and  scattered  electric  fields: 


fihh  =  ES  /  Eih  jES^O 

Shv  =  E^  /  E'v  E'h=0 
£vh  =  Esv  /  E’h  |E5v=0 
£w  =  Esv  /  EV  IeH^O 


(2-8) 


In  general,  eight  quantities  are  needed  to  completely  specify  the 
scattering  matrix.  These  eight  quantities  consist  of  four  magnitudes  and 
four  phases.  The  number  can  be  reduced  to  seven  quantities  if  one  of  the 
phases  is  used  as  a  reference.  If  the  radar  operates  monostatically,  then 
£hv  =  £yh,  and  [  £]  can  be  specified  by  five  quantities  (Knott, 1985).  Once  the 

scattering  matrix  is  determined,  then  all  of  the  scattering  properties  of  the 
target  are  known  for  that  frequency  and  target-radar  orientation. 

For  a  dihedral  corner  reflector  oriented  so  that  its  seam  is  at  an 
angle  a  with  respect  to  the  horizontal  reference  (principal)  plane,  the 
scattering  matrix  for  the  monostatic  case  is  (Cole, 1982) 


[S] 


cos  2a  sin  2a 


1 1 


sin  2a  -cos  2a 


(2-9) 


The  scattered  field  can  then  be  represented  by 


- 1 

IM 

_ 1 

cos  2a 

sin  2a 

E*h‘ 

.  Esv. 

sin  2a 

-cos  2a 

(2-10) 


U..ni  forni-Theoiy .  .of .Diffraction 

The  Uniform  Theory  of  Diffraction  (UTD)  (Kouyoumjian,1974) 
developed  by  Kouyoumjian  and  Pathak  is  valid  for  two  and  three 
dimensional  edge  effects  of  flat  plates  with  straight  edges.  Figure  2.1 
shows  the  geometry  for  the  wedge  diffraction  problem.  A  source  located  at 
CKs'.d'.z’)  radiates  an  incident  electric  field  J^s')  of  arbitrary  polarization, 

which  can  have  a  plane,  cylindrical,  or  spherical  wavefront.  The 
diffraction  of  the  incident  field  occurs  at  the  point  Qe  (or  points  in  the  case 

of  a  plane  wave)  which  is  unique  for  a  given  source  and  observation 
location.  Kouyoumjian  and  Pathak  have  shown  that  the  calculation  of  the 
diffracted  field  at  the  observation  point  P(s,$,z)  can  be  simplified  if  a  fixed 
coordinate  system  centered  at  Qe  is  used  (Kouyoumjian,! 974).  The 


diffraction  of  the  incident  ray  forms  a  cone  of  diffracted  rays  whose  axis  is 
the  edge  with  a  cone  angle  equal  to  the  incident  angle,  i.e.  (J0  =  (5’0. 


i 


t 


12 


1 


The  set  of  orthogonal  unit  vectors  associated  with  the  ray  fixed 
coordinate  system  are  defined  as  follows: 

A  A  A 

I  =  p  o  x 

A  A  A 

D  =  P>0  x  <$) 

A  A 

where  I  is  the  incident  propagation  direction  unit  vector,  and  D  is  the 
diffraction  propagation  direction  unit  vector.  The  diffracted  field  is  given 
by  (Kouyoumjian,1974) 

Ed(s)  -  £i(QE)  •  <Be(D,T)>  A(s)  e#s  (2-11) 

where  the  diffraction  dyadic  is  given  by 

<nE(D,T)>  =  -  P'oBoIIs  -  WDh  (2-12) 

and  the  terms  of  the  diffraction  dyadic  are 

Es.hfL.M'.Po)  =  t  ( -  e'jn/4 )  /  (  2n  (271k)1/2  sinp0  )  ] 
x  { cot  [  (n+p  )  /  (2n)  ]  £  [  kLa+(p-)  ] 

+  cot[(x-p-)/(2n)l  E  [  kLa-(p-)  ] 

*  ( cot  [  (x+p+)  /  (2n)  ]  E  [  kLa+(p+)  ] 

+  cot  [  (JC-P+)  /  (2n)  3  E 1  kLa-(P+)  ] ) }  (2-13) 
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where 


L  =  distance  parameter  dependent  on  wavefront 
k  =  wavenumber  (2 %  /  wavelength  ) 

P±  =  sum/difference  angle  of  D  and  T  as  measured  in  the  x-y  plane 
P0  =  angle  of  reflection  as  measured  from  the  z  axis 
A(s)  =  spread  factor  of  wavefront 
n  =  wedge  angle  number 

and 

a±(P)  =  2  cos2  [  { 2nxN±  -  (p) )  /  2  ]  (2-14) 

where  N±  are  the  integers  which  most  nearly  satisfy  the  equations 
2nxN+  -  (P)  =  7i 

and 


2mtN*  -  (p)  =  -x 


with 


P  =  p1  =  $  ±  <{>’ 

E  is  a  transition  function  that  was  included  by  Kouyoumjian  and 
Pathak  to  improve  the  diffracted  field  solution  in  the  transition  regions  of 
the  incident  and  reflection  shadow  boundaries.  E  (x)  is  defined  by 


E(x)  -  2j  |x1/2j  ei*  j  e*Jt2  dt 


(2-15) 


I 


The  magnitude  and  phase  of  £(x)  are  plotted  in  Figure  2-2.  Kouyoumjian 
and  Pathak  used  the  following  approximation  for  small  arguments  (x<0.3) 

E(x)  ~  [Ux)l/2-2xej*/M2/3)x2e-jK/4]ej(*/4  +  *)  (2-16) 

and  for  large  arguments  (x>5.5) 

E(x)~  1  +  j  /  ( 2  x )  -  3  /  (4x2 )  -  j  15  /  ( 8  x3 )  +  75  /  (16X4)  (2-17) 

If  n=2,  then  the  wedge  reduces  to  a  half  plane  and  the  diffraction 
coefficient  reduces  to 

Ds.ha.M'.po)  =  [  ( -  e  -i  "/< )  /  (  2  {2xk)i/2  sinp0 )  ] 
x[E(kLa(p-)l/cos(  P-/2) 

*  E  { kLa  (p+) )  /  cos(  P+  /  2  )  ]  (2-18) 

Another  special  case  of  the  diffraction  coefficients  occurs  when  the 
source  and  receiver  are  in  the  far  field  of  the  wedge,  so  that  E00  is 
approximately  equal  to  1 .  For  this  case  the  diffraction  coefficients  reduce  to 
Keller's  diffraction  coefficients  (Keller,*  962): 

Q,,h(L,M\p0)  =  [  ( -  e  -  J  */< )  /  ( 2  {2xk)l/2  sin[50 )  ] 

x  [  cot(  {  n  +  P* )  /  (2n) )  +  cot(  { x  -  p  1  /  {2n} ) 

*  { cot(  { x  +  P+ }  /  {2n} )  +  cot(  { x  -  P+ }  /  {2nl ) )  ]  (2-19) 
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FIGURE  2-2.  TRANSITION  FUNCTION 
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The  diffracted  fields  can  be  written  in  matrix  form  as 

Edp0(s)  _Els  0  E^Qe) 

=  A(s)e-iks  (2-20) 

.  Ed*(s)  J  L  0  -HhJ  LEV(Qe) 

The  diffraction  dyadic  matrix  is  diagonal  because  of  the  ray  fixed 
coordinate  system  employed.  The  subscripts  s  and  h  refer  to  the  soft  and 
hard  boundary  conditions  respectively.  The  soft  diffraction  coefficient,  Hs> 

uses  the  soft  boundary  condition  (E  fieid  component  parallel  to  the  edge): 

Ell  I  wedge  =  0  (2-21) 

and  the  hard  diffraction  coefficient,  Eh,  uses  the  hard  boundary  condition 
(E  field  component  perpendicular  to  the  edge): 

9E_l  /  ^  j  wedge  =  0  (2-22) 
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The  spread  factor,  A(s),  and  the  distance  parameter,  L,  are 


dependent  on  the  type  of  incident  wave  such  that 


A(s)  =  \ 


(S)  -1/2 


(S*  /  (S  (s+s'))1/2 


plane,  cylindrical,  &  conical 
wave  incidence  (2-23) 

spherical  wave  incidence 


L  = 


s  sin2  p0 
•*  s  s'  /  (s+s') 
jss'  (sin2  P0)  /  (s+s') 


plane  wave  incidence 

cylindrical  wave  incidence  (2-24) 

conical  &  spherical  wave  incidence 


The  Method  of  Moments 

The  Method  of  Moments  (MOM)  is  a  general  procedure  for  solving 
linear  equations  (Harrington, 1 966,1 968).  The  MOM  starts  with  a  general, 
inhomogeneous  equation: 

L(f)  =  g  (2-25) 

where  L  is  an  arbitrary,  linear  operator,  g  is  the  known,  and  f  is  the 

unknown  to  solve.  Begin  by  letting  f  be  represented  by  a  series  solution  of 
fn,  where  fn  is  in  the  domain  of  the  operator  L: 


f  =  Zanfn  (2-26) 

where  an  are  the  constant  coefficients  to  find.  The  functions  fn  are  called 
the  expansion  or  basis  functions.  Usually  an  infinite  summation  of  a 
complete  set  of  orthogonal  expansion  functions  is  needed  for  an  exact 
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solution,  but  only  a  finite  summation  is  needed  for  an  approximate 
solution. 


Now  substitute  Equation  (2-26)  into  (2-25)  to  obtain 

L  ( I  a„  f„  )  =  g  (2-27) 

and  use  the  linearity  property  of  L  to  get 

XanL(fn)  =  g  (2-28) 

Next,  introduce  an  inner  product  <x,y>,  which  is  a  scalar  defined  by 

<x,y>  =  (  x(t)  y(t)  dt  (2-29) 

0 

A  set  of  weighting  or  testing  functions  wm  must  be  defined  in  the  range  of 
L.  Take  the  inner  product  of  Equation  (2-28)  with  w^  to  get 


T  an  <wm,L(fn)>  =  <wm,g>  (2-30) 

With  m  =  1,2,3,***,  a  system  of  equations  can  now  be  written  in  matrix  form: 
t  U  U  a„  ]  =  [  gm  ]  (2-31) 

where 


[1 


mn 


1  = 


<w1,L(f1)>  <wltL(f2)>  •  •  • 

<w2,L(fi)>  <w2,L(f2)>  *  *  * 


(2-32) 


L 
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<W1(g> 

<w2,g> 

[gm]=  •  (2-33) 

As  long  as  the  inverse  of  [  lmn  ]  exists,  then  an  is  given  by 

[«„]  -  [lnmHUml  <2*34> 

Representing  fn  as  a  matrix  of  functions 

[  fn  3  =  [  fi  f2  f3  *  *  ]  (2-35) 

the  solution  of  f  may  be  written  as 

f  =  [f„][a„]  =  tf„Hlnm]1[gm]  (2-36) 

This  solution  may  be  approximate  or  exact,  depending  on  the  choice 
of  the  expansion  (basis)  and  testing  (weight)  functions,  fn  and  wm 

respectively.  A  special  case  called  the  Galerkin’s  Method  occurs  when 
fn=wn  • 
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Next,  applying  the  MOM  to  an  electromagnetic  field  problem,  start 
with  Maxwell's  equations  in  the  time  harmonic  form: 

Vx£  =  -jwpH 

VxH  =  jtoeE  +  JL  (2-37) 

where  E  is  the  electric  field  intensity,  H  the  magnetic  field  intensity,  and  jl 

the  electric  current  density.  Combining  the  above  equations  into  a  single 
equation  for  E  yields 

-(jco)-i  Vx(p-iV  xE)-jcoeE  =  si  (2-38) 

In  the  form  of  the  MOM  operator  notation,  this  becomes 

-L(E)  =  iI  (2-39) 

Consider  an  arbitrarily  shaped  scatterer  in  a  homogeneous 
medium.  An  impressed  source  (  iL‘ ,  M.’  )  radiates  the  fields  (  E' .  11* )  in 
free  space  and  the  fields  (  ET  t  H?  )  in  the  presence  of  the  scatterer.  For  a 
perfect  electric  conductor,  surface  currents  are  induced  on  the  surface 

of  the  scatterer,  which  produce  the  scattered  field  Es-  Apply  Equation  (2- 
39)  to  the  scatterer  to  get 


-L(ES)  =  iZs 


(2-40) 


Applying  the  inverse  operator  to  both  sides  yields  the  potential 
integral  solution  to  Equation  (2-40): 

-Es  =  L'x(ils)  =  -jo>A8-V$s  (241) 

where 

As(r)  =  |i  ||{  jls(r')  (4rclr- r'l)-1  e-ik,r-r’1  dv' 

3>s(r)  =  e*1  {{ j  Cs(r')  (4ji  I  r  -  r’  I  )-1  e'ik  1  r  • r’ 1  dv' 

as  =  *  ( jco  )-1  v  •  jig 

On  the  surface  of  the  scatterer,  the  boundary  conditions  on  the 
electric  field  are  such  that  the  tangential  component  must  vanish: 


nx£T  =  n  x  ( E’  +  Es )  =  0  (242) 

Thus 

Eitan  =  *Estan  (2-43) 

where  the  subscript  "tan"  means  the  tangential  component  of  the  electric 
field.  Thus,  Equation  (241)  becomes 

E'tan  =  L'1  (  )  (244) 
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To  continue  along  the  lines  of  the  MOM  solution,  a  suitable  inner 
product  must  be  chosen.  For  electromagnetic  field  problems,  a  suitable 
inner  product  is 


<£,*!> 


E  •  *1  dv' 


(2-45) 


which  is  the  definition  of  reaction  (Harrington, 1 966,1 968).  Equation  (2-45) 
can  be  expressed  as  an  integral  over  the  surface  of  the  scatterer  since  jl  is 
the  surface  current  density  and  E  is  limited  to  the  tangential  component  on 
the  surface. 

The  next  step  is  to  expand  the  surface  current  density  ils  in  terms  of 
N  expansion  (basis)  functions  tin- 


i  =  Ilniin  (246) 

Substitute  Equation  (246)  into  Equation  (2-44)  to  get 

L'MIlnJn)  »EU  (247) 

and  since  L  is  a  linear  operator,  Equation  (2-47)  becomes 

I  In  L-1  <dn>  =  E’tan  (2-48) 

Next  choose  a  set  of  testing  (weight)  functions  2Ln  which  are  also 

tangentir.I  10  the  surface.  Apply  the  inner  product  defined  in  Equation  (2- 
45)  to  Equation  (2-48)  and  to  obtain 

I  In  (*,)>  =  <3Em£it.„>  (249) 
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Now  the  system  of  equations  can  be  represented  in  matrix  notation: 


[  Zmn  ] 

[Ini 

=  [Ym] 

where 

II 

OEiE) 

12 

<5£2.E'> 

[Inl  = 

[Ym]  = 

• 

• 

- 

<SiX‘l 

il2)>  •  •  • 

[  2mn  ] 

<5E2,L1  (Ji)> 

<S2,L-1(J2)>  •  .  . 

Solving  for  [ 

]  one 

obtains 

tin]  = 

Z„m  J1  [  Ym  ] 

and 


>1  =  =  [41J[Z„,n]-1m,] 


f 

r 


(2-50) 


(2-51) 


(2-52) 


(2-53) 
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where  the  elements  of  the  impedance  and  excitation  matrix  are  given  by 
^mn  =  "  \\  Em  '  sin 

Yn.  -  Ui  in,  •  £'  dv  (2-54) 

Zmn  is  integrated  over  the  surface  of  the  n-th  expansion  mode,  and  is 
called  the  mutual  impedance  between  modes  m  and  n.  is  integrated 

over  the  volume  occupied  by  the  source  (  il’ ,  M.’  )»  and  is  called  the  modal 
excitation  voltage. 
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III.  Analysis 


Ch<?pging  the.,  Approach 

During  the  initial  discussions  with  Zelnio,  the  sponsor  of  this 
research,  the  approach  to  attacking  the  problem  of  the  dihedral  comer 
reflector  could  have  taken  several  paths.  The  sponsor’s  contract  for  the 
indoor  measurements  was  with  The  Ohio  State  University  ElectroScience 
Laboratory  (OSU-ESL)  in  Columbus  Ohio,  so  one  obvious  approach  was  to 
use  the  Uniform  Geometrical  Theory  of  Diffraction  (UTD),  which  was 
greatly  advanced  by  OSU.  For  comparison  and  validation  purposes, 
another  approach  should  be  pursued  to  check  against  measurements. 

When  used  properly,  the  Method  of  Moments  (MOM)  approaches  the 
exact  solution  of  Maxwell's  integral  equations.  Therefore,  the  MOM  was 
thought  of  as  the  best  approach  to  compare  with  the  UTD  results  and  the 
measurements.  The  problem  with  MOM  is  that  as  the  object  modeled 
becomes  larger  in  electrical  size  (size  in  wavelengths)  the  amount  of 
computer  memory  required  increases,  and  the  CPU  time  increases 
substantially.  Thus,  the  computer  resources  limit  the  size  of  an  object 
which  could  be  modeled  by  the  MOM.  Much  research  has  been  done  on 
modifying  the  MOM  to  handle  larger  problems.  From  this  research 
several  special  techniques  have  emerged  which  could  solve  MOM  problems 
faster  and  handle  larger  electrical  sizes. 

Since  a  9  inch  dihedral  comer  reflector  at  10  gigahertz  (GHz)  is  116 
square  wavelengths  in  electrical  size,  an  alternative  to  the  original  MOM 
approach  should  be  investigated.  Even  at  3  GHz  the  dihedral  is  10.5  square 
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wavelengths,  which  would  be  considered  a  large  problem  for  MOM. 
Should  one  of  the  special  MOM  techniques  be  pursued  or  should  an 
alternative  solution  be  pursued? 

After  a  thorough  search  of  the  literature,  one  technique  seemed  to  be 
the  most  promising  for  the  problem  at  hand.  The  Spectral  Domain 
Approach,  also  called  the  Spectral  Theory  of  Diffraction  (STD) 
(Mittra,1972),  was  developed  by  Mittra  and  Li  at  the  University  of  Illinois  at 
Urbana-Champaign.  The  STD  uses  a  UTD  solution  as  an  initial  guess, 
and  then  through  an  iterative  technique  comes  to  a  final  solution. 

The  Southeastern  Center  for  Electrical  Engineering  Education 
(SCEEE)  sponsored  a  short  course  on  RCS  computational  techniques  at  the 
University  of  Illinois  at  Urbana-Champaign,  and  Mittra  was  the 
coordinator  of  the  short  course.  At  the  short  course,  the  problem  of  the 
dihedral  corner  reflector  was  presented  to  Mittra.  The  STD  has  been  used 
to  obtain  results  on  a  square  flat  plate  and  on  an  infinite  rectangular 
cylinder,  both  of  which  are  related  to  the  dihedral  comer  reflector  problem. 
However,  Mittra  suggested  that  the  STD  would  be  too  difficult  to  apply  to 
the  dihedral  problem,  and  that  there  were  other  techniques  which  would 
be  better  suited  for  the  dihedral. 

Mittra  proposed  using  a  special  technique  with  the  MOM.  His 
suggestion  involved  using  the  Physical  Optics  (PO)  approximation  of  the 
surface  current  as  an  entire  domain  expansion  (basis)  function  on  the 
portions  of  the  dihedral  surface  which  were  sufficiently  far  enough  away 
from  the  edges.  On  the  regions  near  the  edges,  the  technique  of 
subsectional  basis  functions  could  then  be  used  with  the  MOM. 
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After  the  investigation  of  the  methods  available  to  attack  the  dihedral 
comer  reflector  problem,  another  factor  had  to  be  considered-time. 
Because  of  the  limited  amount  of  time  alloted  to  complete  the  research,  an 
approach  or  approaches  had  to  be  selected  on  the  merit  that  they  could  be 
completed.  Aberegg  and  Marhefka  at  OSU  had  been  working  on  a 
computer  code  utilizing  the  UTD  to  solve  a  two  dimensional  dihedral 
corner  reflector.  This  seemed  to  be  the  logical  choice  for  doing  the  UTD 
analysis  on  the  comer  reflector.  Next,  the  alternative  method  needed  to  be 
chosen. 

With  the  knowledge  that  a  Cray  supercomputer  was  available  for 
computation,  a  general  MOM  solution  was  sought  after  which  could  be 
modified  for  the  dihedral.  The  supercomputer  would  alleviate  the  problem 
of  the  slowness  of  the  MOM  solution.  A  general  purpose  MOM  code,  the 
Electromagnetic  Surface  Patch  (ESP)  Code:  Version  II  -  Polygonal  Plates 
and  Wires  (Newman, 1985),  was  found  through  the  help  of  the  people  at 
OSU.  The  ESP  code  is  based  on  the  Electric  Field  Integral  Equation  (EFIE), 
so  it  was  ideal  for  open  structures  such  as  the  dihedral  comer  reflector. 

Now  that  a  UTD  code  and  a  MOM  code  had  been  obtained,  the  attack 
on  the  dihedral  corner  reflector  problem  could  begin.  The  plan  was  to  run 
the  UTD  corner  reflector  code  on  a  VAX  11/785  computer,  and  since  the 
code  was  developed  (and  still  being  developed  at  the  time  of  this  writing)  on 
a  VAX  11/780  computer  this  should  have  been  no  problem.  The  ESP  code 
was  also  developed  on  the  VAX  11/780  computer,  but  the  intention  of 
running  the  code  on  the  Cray  X-MP  12  Serial  315/75  computer  opened  up 
the  potential  trouble  area  of  Fortran  compiler  ambiguities. 
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Because  of  the  potential  trouble  area  with  the  Fortran  compilers,  the 
implementation  of  the  ESP  code  on  the  Cray  was  the  first  approach  tried  on 
the  dihedral  comer  reflector.  The  Cray  needs  another  computer  connected 
to  it  as  a  front  end  for  obtaining  interactive  access.  This  front  end 
computer  was  also  a  VAX  11/785.  After  the  ESP  code  was  uploaded  from 
the  front  end  VAX,  the  testing  of  the  ESP  code  on  the  Cray  began. 

The  first  thing  the  Cray  didn't  recognize  was  the  use  of  special 
control  characters  in  the  code  which  performed  various  functions  such  as 
tabbing,  line  feeding,  and  form  feeding.  Some  of  the  control  characters 
were  removed  globally  with  the  VAX  editor,  while  others  had  to  be  removed 
on  a  line  by  line  basis.  Considering  the  fact  that  the  ESP  code  was  over  7800 
lines  in  length,  stripping  all  of  the  special  control  characters  was  no  quick 
and  easy  task.  Once  the  special  control  characters  were  removed,  more 
Fortran  compiler  ambiguities  became  apparent.  The  most  frequently 
encountered  difference  was  in  the  use  of  input/output  (I/O)  statements, 
such  as  "open"  and  "close"  statements.  After  modifying  the  code  to  the 
format  required  by  the  Cray,  it  finally  compiled  without  errors. 

The  code  operated  with  an  input  parameter  file  which  needed  to  s 
assigned  to  the  logical  device  unit  number  "5".  This  unit  number  was  then 
opened  by  the  ESP  code  and  the  program  read  the  input  parameter  file. 
After  some  final  fine  tuning,  the  ESP  code  was  ready  to  be  tested  on  the 
Cray.  Starting  with  input  parameters  provided  by  Newman,  the  author  of 
the  ESP  code,  the  version  of  the  code  on  the  Cray  could  be  validated  with 
known  results.  The  Cray  performed  superbly  in  calculating  the  test  run 
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approximately  20  times  faster  than  the  VAX.  How  well  would  the  Cray 
perform  on  the  dihedral  comer  reflector  at  3  GHz? 

As  stated  previously,  the  dihedral  at  3  GHz  is  approximately  10.5 
square  wavelengths  electrically.  Sampling  4  times  per  wavelength  yielded 
approximately  24  modes  per  square  wavelength.  Including  overlap  modes, 
this  gave  a  total  of  261  modes  for  which  to  solve.  The  ESP  code  was 
originally  dimensioned  to  allow  up  to  500  total  modes,  so  half  of  the  total 
was  used.  For  better  accuracy,  sampling  around  8  times  per  wavelength 
resulted  in  112  modes  per  square  wavelength  for  a  total  of  1194  modes, 
which  was  more  than  double  the  dimensions  in  the  ESP  code,  so  the 
dimensions  of  the  arrays  in  the  code  needed  increasing. 

As  a  first  attempt,  the  arrays  were  increased  to  2000,  four  times  the 
original  value.  The  code  compiled  without  problems,  but  during  execution 
the  Cray  aborted  the  run  with  an  indication  that  more  memory  was 
requested  than  available.  After  consultation  with  the  Cray  Research 
analysts,  the  conclusion  was  drawn  that  the  Cray  couldn’t  handle  the  code 
with  the  arrays  dimensioned  that  large.  The  Cray  X-MP  12  Serial  315/75 
has  16  megabytes  (MB)  of  main  memory,  but  each  word  is  8  bytes  long  so 
the  machine  has  a  total  of  2  megc  vords  (MW)  of  main  memory.  The  code 
itself  used  about  30  kilowords  (KW),  and  the  Cray  Operating  System  (COS) 
consumed  about  200  KW  of  memory,  leaving  1 .77  MW  for  the  variables  in 
the  code. 

After  careful  examination  of  the  code  and  consultation  with 
Newman,  a  dimension  of  about  1000  appeared  to  be  the  largest  the  Cray 
could  handle.  With  the  code  dimensioned  to  allow  arrays  of  up  to  1000 
elements,  the  code  compiled  and  executed  without  trouble.  Thus,  the 
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maximum  the  dihedral  could  be  sampled  at  3  GHz  was  7  times  per 
wavelength. 

At  10  GHz  and  with  the  arrays  dimensioned  at  1000,  the  maximum 
the  dihedral  could  be  sampled  would  have  been  only  two  times  per 
wavelength!  Not  nearly  enough  even  for  an  approximate  result,  and  an 
accurate  answer  was  desired.  At  this  point,  it  became  very  apparent  that 
the  Cray  was  not  the  solution  to  running  a  large  MOM  program.  The  Cray 
has  the  computational  speed  required,  but  lacks  the  required  main 
memory  for  a  large  matrix  problem.  The  Cray  analysts  were  aware  of  the 
memory  problem  and  had  plans  of  upgrading  the  memory  to  48  MB  (3 
times  the  current  memory)  in  May  of  1987.  However,  this  planned  memory 
increase  didn't  help  the  current  problem,  and  would  not  be  implemented 
in  time  for  this  thesis. 

Now  there  was  a  big  problem-how  to  run  the  ESP  code  for  larger 
objects.  Turning  to  the  only  other  computer  at  hand,  the  VAX  11/785,  the 
version  of  the  ESP  code  with  the  2000  dimensioned  arrays  was  compiled, 
but  couldn't  be  linked  because  of  a  lack  of  virtual  address  space.  The  VAX 
system  operator  increased  the  virtual  address  space,  but  once  again  the 
program  wouldn't  link  due  to  the  pagefile  quota  being  exceeded.  After  the 
pagefile  was  increased  by  the  system  operator,  the  program  linked  without 
any  trouble,  but  upon  execution  of  the  program  the  VAX  aborted  the  run 
due  to  the  memory  being  full.  The  VAX  system  message  and  error 
recovery  documentation  was  very  useful,  and  led  to  a  solution  to  increase 
the  memory  that  any  single  job  could  use.  This  was  defined  by  the  working 
set  default  (WSdef)  and  extent  (WSext)  in  both  the  user  authorization  file 
(UAF)  and  the  system  generation  utility  (SYSGEN)  With  all  the  changes 
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made  to  the  operating  system,  the  program  finally  executed  for  the 
dihedral  at  3  GHz  with  7  samples  per  wavelength. 

With  all  the  work  required  to  run  the  program  at  3  GHz,  what  would 
it  take  to  get  the  program  to  run  at  10  GHz?  Again,  the  virtual  address 
space,  page  file  quota,  and  the  working  set  parameters  were  increased. 
The  dimensions  in  the  ESP  code  were  now  increased  to  4000,  allowing  the 
dihedral  to  be  analyzed  at  10  GHz  with  a  sampling  rate  of  4  times  per 
wavelength,  the  minimum  needed  for  engineering  accuracy.  Again,  the 
VAX  aborted  the  job  because  of  a  lack  of  memory.  This  time  it  was  thought 
that  the  Digital  Electronic  Corporation  (DEC),  the  makers  of  the  VAX, 
should  be  called  in  on  this  problem. 

The  DEC  analyst  for  the  VAX  suggested  configuring  the  system  to 
allow  the  user  full  access  of  the  computer's  resources.  Using  a  backup  of 
the  operating  system,  the  system  disk  was  configured  to  the  specifications 
the  VAX  analyst  suggested.  This  time  the  program  began  to  execute.  To 
the  casual  observer,  it  appeared  that  the  VAX  was  operating  normally. 
However,  the  next  morning  the  system  slowed  down  to  a  standstill,  and  the 
system  had  to  be  powered  down  and  rebooted.  Since  the  trouble  didn't 
appear  to  begin  until  other  users  were  on  the  system,  and  knowing  that  the 
program  was  allowed  the  entire  system's  resources,  it  seemed  possible 
that  the  program  might  run  over  the  weekend  when  there  were  no  other 
users  with  which  to  compete. 

With  the  anticipation  that  the  program  would  execute,  it  was  again 
submitted.  After  about  18  hours  of  CPU  time,  it  became  apparent  that  the 
system  had  slowed  down  to  a  standstill  once  again.  But  a  message  on  the 
operator's  console  provided  insight  as  to  why  the  program  wouldn't  run. 
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Since  the  VAX  is  a  virtual  addressing  computer,  it  pages  in  and  out  of  the 
main  memory  to  a  disk  file  called  the  pagefile.  The  problem  was  that  the 
pagefile  was  not  large  enough  to  handle  all  the  pages  created  by  the 
program.  This  problem  was  alleviated  by  increasing  the  size  of  the 
pagefile  to  approximately  four  times  the  size  of  the  original  pagefile.  This 
allowed  the  program  space  to  page  in  and  out  of  the  main  memory.  The 
ESP  code  finally  started  executing  without  signs  of  trouble.  After  12  days 
elapsed,  the  VAX  had  to  be  powered  down  due  to  a  problem  with  the  air 
conditioning  system.  At  this  point,  the  code  had  used  approximately  120 
hours  of  CPU  time,  and  had  completed  77%  of  the  360  total  field  points. 
Because  of  the  dihedral's  symmetry  in  the  monostatic  case,  the  remaining 
points  could  be  extrapolated. 


Implementation  of  the  UTD  Code 

When  it  became  apparent  that  there  would  be  problems  running  the 
ESP  code,  the  analysis  of  the  dihedral  comer  reflector  was  started  with  the 
UTD  two  dimensional  comer  reflector  code  developed  by  Aberegg  and 
Marhefka  at  OSU.  The  code  included  first,  second,  and  third  order  effects, 


which  included  single,  and  double  diffract  n.  Thus,  the  code  included  the 
following  terms: 

single  diffracted  (D) 

diffracted-reflected  and  reflected-diffracted  (DR,RD) 
reflected-diffracted-reflected  (RDR) 

diffracted-double  reflected  and  double  reflected-diffracted  (DRR,  RRD) 
double  diffracted  (DD) 
diffracted-reflected-diffracted  (DRD) 

double  diffracted-reflected  and  reflected-double  diffracted  (DDR,  RDD) 
The  code  was  a  two  dimensional  model,  so  it  couldn’t  predict  cross 
polarization  components  or  off  principal  plane  scattering. 

A  very  strong  advantage  of  the  UTD  solution  is  the  speed  of 
computation  and  the  small  amount  of  computer  resources  required.  Since 
the  diffraction  coefficients  consist  of  little  more  than  a  combination  of 
trigonometric  functions,  solutions  are  easily  and  speedily  obtained.  This 
allowed  the  dihedral  to  be  analyzed  at  10  GHz,  where  the  ESP  code  had  a 
difficult  time  doing  because  of  computer  limitations.  Another  feature  of 
the  UTD  analysis  is  the  RCS  is  broken  up  into  the  individual  components  of 
the  scattering  interactions,  thus  bringing  insight  into  the  physical 
phenomena. 

Unlike  the  ESP  code,  the  UTD  code  ran  interactively,  with  the  user 
entering  responses  after  the  code  prompted  for  an  input.  The  program 
first  asks  for  the  width  of  the  strips  and  the  flare  angle  between  them. 
Since  the  code  was  a  two  dimensional  model,  the  dihedral  comer  reflector 
only  had  width  as  a  variable,  with  the  length  being  infinite.  The  code  then 
prompts  for  either  a  monostatic  or  bistatic  case.  If  the  run  is  for  a  bistatic 
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case,  then  it  prompts  for  the  bistatic  angle.  The  next  7  prompts  are  for  the 
scattering  terms  the  code  is  to  include  in  the  run.  Lastly,  the  code  prompts 
for  either  the  soft  or  hard  boundary  condition.  Even  if  all  the  scattering 
effects  were  included,  the  code  only  took  a  few  seconds  to  complete,  unlike 
the  MOM  code  which  took  hours. 

The  UTD  comer  reflector  code  worked  well.  The  only  problem  areas 
occurred  when  the  observation  point  was  on  the  diffraction  shadow 
boundary  (SB)  or  reflection  shadow  boundary(RB).  The  DD  component 
tended  toward  infinity  rather  than  approaching  the  proper  limit.  A 
uniform  double  diffraction  solution  had  only  recently  been  solved  by  Tiberio 
and  Kouyoumjian  (Aberegg,1985)  and  was  not  yet  implemented  in  this 
code. 

The  UTD  corner  reflector  code  did  not  experience  any  problems  of 
the  magnitude  that  the  ESP  code  experienced.  The  ESP  code  presented 
quite  a  demand  for  computer  system  resources,  and  as  a  result  required 
more  time  and  effort  in  running.  However,  the  ESP  solved  for  both  cross 
and  same  sense  polarization,  so  it  was  preferred  when  the  object  was  not 
electrically  large. 
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In  this  chapter  the  results  of  the  UTD  two  dimensional  corner 
reflector  code  and  the  ESP  code  will  be  presented  and  compared  to 
measurements  taken  at  The  Ohio  State  University  ElectroScience 
Laboratory  (OSU-ESL).  The  geometry  of  the  dihedral  corner  reflector  is 
shown  in  Figure  4-1.  The  dihedral  consists  of  two  9  by  9  inch  square  flat 
plates  joined  to  form  an  interior  angle  of  90  degrees.  The  plates'  material 
is  assumed  to  be  perfectly  conducting. 
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Figure  4.1  Dihedral  Corner  Reflector 

The  dihedral  comer  reflector  was  analyzed  at  2,  3,  and  10  GHz  using 
the  UTD  two  dimensional  comer  reflector  code.  Because  the  UTD  is  a  high 
frequency  approximation,  the  accuracy  of  the  results  is  expected  to  degrade 


as  the  the  frequency  is  decreased.  The  ESP  moment  method  code  was  used 
to  analyze  the  dihedral  at  2  and  3  GHz  with  the  possibility  of  obtaining 
results  at  10  GHz.  The  accuracy  of  the  ESP  code  does  not  depend  on 
frequency,  but  rather  on  the  wavelength  sampling  rate.  However,  as  the 
frequency  is  increased  the  matrix  in  the  ESP  code  becomes  large  which 
forces  the  sampling  rate  to  be  lowered,  thus  affecting  the  accuracy. 

The  dihedral  was  modeled  for  both  monostatic  and  90°  bistatic 
geometries.  The  dihedral  was  also  modeled  off  the  principal  plane  with  the 
ESP  code.  Recall  that  the  UTD  cannot  model  the  dihedral  off  the  principal 
plane  because  the  model  is  two  dimensional.  Cross  polarization  could  only 
be  modeled  with  the  ESP  code,  because  the  UTD  code  is  two  dimensional. 

ESP  Code  Solutions 

At  what  wavelength  sampling  rate  does  the  ESP  code’s  solution 
converge  when  analyzing  the  dihedral?  What  testing  function  should  be 
used  to  achieve  the  most  accurate  results?  What  are  the  tradeoffs  between 
the  two  testing  functions  in  terms  of  accuracy,  computation  speed,  and 
computer  storage?  These  questions  need  to  be  considered  when  using  the 
ESP  code  to  model  the  scattering  from  any  target. 

The  solution  the  ESP  code  generates  becomes  more  accurate  as  the 
sampling  rate  is  increased  per  wavelength.  Figure  4-2  shows  the 
monostatic  RCS  of  the  dihedral  at  3  GHz  for  theta  in  theta  out  polarization 
(soil  boundary  condition  on  the  electric  field),  and  for  phi  in  phi  out 
polarization  (hard  boundary  condition  on  the  electric  field).  The  graph 
shows  the  comparison  between  filament  testing  and  full  surface  patch 
testing  with  a  0.25  wavelength  patch  (4  samples  per  wavelength).  For  theta 
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polarization,  the  two  methods  differ  slightly  in  the  0°  to  20°,  70°  to  130°,  and 
150°  to  180°  regions.  The  most  significant  discrepancy  occurs  around  110° 
where  the  location  of  the  nulls  differ.  Despite  the  minor  differences  the  two 
methods  show  excellent  agreement. 

Figure  4-3  compares  the  two  testing  methods  with  a  0.20  wavelength 
patch  (5  samples  per  wavelength).  There  is  only  slight  disagreement 
between  the  two  methods  in  the  same  aspect  angle  regions,  but  the 
differences  are  decreasing.  In  Figure  4-4  the  two  methods  again  show 
slight  disagreement,  but  are  converging  to  the  same  result  with  patch  size 
of  0.15  wavelengths  (6.67  samples  per  wavelength).  The  most  significant 
difference  occurs  in  the  90°  to  100°  region  of  the  phi  in  phi  out  pattern,  but 
differs  by  2  dB  at  most.  This  gives  reasonable  confidence  in  the  two 
separate  testing  methods  as  they  both  tend  toward  the  same  result. 

Figure  4-5  shows  the  comparison  between  0.15,  0.20,  and  0.25 
wavelength  patches  for  filament  testing,  and  the  largest  difference  is  at 
110°  for  theta  in  theta  out.  Elsewhere  there  is  good  agreement  between  the 
sampling  rates  for  filament  testing.  Figure  4-6  shows  the  same 
comparison  for  full  surface  testing,  For  theta  in  theta  out,  the  worst 
agreement  occurs  at  110°,  a  difference  of  almost  5  dB.  But  elsewhere  only 
small  differences  exist  between  the  sampling  rates  giving  confidence  in  the 
results,  with  the  0.15  wavelength  patch  giving  the  most  confidence. 
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FIGURE  4-3.  MONOSTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 

FILAMENT  VERSUS  FULL  SURFACE  TESTING  COMPARISON 
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FIGURE  4-4.  MONOSTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 

FILAMENT  VERSUS  FULL  SURFACE  TESTING  COMPARISON 
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FIGURE  4-5.  MONOSTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 

PATCH  SIZE  COMPARISON  FOR  FILAMENT  TESTING 


Noting  that  full  surface  testing  yields  the  most  accurate  solution,  a 
worst  case  comparison  is  made  in  Figure  4-7.  Filament  testing  is  used 
with  a  0.25  wavelength  patch,  and  full  surface  testing  is  used  with  a  0.15 
wavelength  patch.  For  theta  in  theta  out,  the  greatest  difference  of 
approximately  5  dB  is  at  110°,  and  for  phi  in  phi  out,  the  greatest  difference 
of  about  5  dB  is  at  89°.  There  is  good  agreement  elsewhere,  especially 
considering  the  fact  that  this  is  a  worst  case  comparison. 

UTD  Code  Solutions 

The  UTD  code  used  is  a  two  dimensional  model,  and  solves  up  to 
third  order  effects  including  at  most  double  diffraction.  All  of  the 
diffraction  and  reflection  mechanisms  up  to  the  third  order  are  included 
when  either  the  expression  'DD  Term  Included’  or  'UTD  with  All  Terms’ 
is  used,  and  all  except  the  double  diffraction  mechanism  are  included 
when  ’DD  Term  Excluded’  is  specified.  When  the  interior  angle  of  the 
dihedral  is  90°,  the  code  does  not  include  diffraction  at  the  junction  of  the 
two  plates,  thus  the  only  diffraction  effects  are  from  the  two  outer  edges  of 
the  dihedral  as  shown  in  Figure  4-8. 


Figure  4-8  Diffraction  Mechanisms 
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FIGURE  4-7.  MONOSTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 

FILAMENT  VERSUS  FULL  SURFACE  TESTING  COMPARISON 
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Figure  4-9  shows  a  double  diffraction  term  comparison  at  2  GHz  for 
both  theta  in  theta  out  and  phi  in  phi  out  polarization  for  the  monostatic 
case.  At  135°  and  180°  for  theta  in  theta  out,  the  solution  diverges  to  infinity 
when  all  the  terms  are  included.  The  solution  blows  up  at  135°  because  the 
observation  point  lies  on  the  reflection  shadow  boundary  (RB)  of  one  of  the 
plates  as  seen  in  Figure  4-10.  Essentially,  a  flat  plate  is  seen,  and  if  the 
code  could  realize  that  fact,  then  the  solution  would  converge.  The  problem 
is  the  code  is  still  using  the  diffraction  terms  from  the  outer  edges  of  the 
two  plates. 
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Figure  4-10  Dihedral  Reflection  Boundary 


When  the  observation  point  approaches  a  reflection  boundary  (RB), 
the  diffraction  coefficients  tend  toward  infinity.  When  a  flat  plate  is  seen, 
the  diffraction  at  the  two  edges  have  opposite  signs,  thus  the  infinities 
cancel  and  the  solution  converges.  But  the  code  sees  one  edge  at  the  RB 


and  the  edge  of  the  second  plate  is  seen  edge  on,  thus  giving  one  infinite 
and  one  finite  coefficient.  To  overcome  this  problem,  the  code  should  treat 
the  diffraction  at  the  junction  of  the  plates  as  a  second  edge  of  the  plate 
with  the  infinite  diffraction  coeffi  This  would  insure  convergence  of 

the  solution. 

At  180°  the  solution  diverges  for  a  similar  reason.  The  problem  at 
180°  is  with  the  double  diffraction  term.  The  second  edge  lies  on  the  RB  of 
the  first  edge,  and  the  observation  point  lies  on  the  RB  of  the  second  edge  as 
seen  in  Figure  4-11. 


I 

x 

Figure  4-11  Dihedral  DD  Reflection  Boundary 


This  problem  is  not  as  easy  to  alleviate  as  the  reflection  boundary 
problem  at  135°.  A  uniform  double  diffraction  coefficient  that  is  well 
behaved  across  the  reflection  boundary  is  needed.  Without  the  DD  term  the 
solution  converges  at  180°  as  seen  in  Figure  4-9.  However,  excluding  the 


DD  term  is  not  without  consequence.  Accuracy  is  lost  throughout  the 
pattern,  and  not  just  in  the  region  near  180°. 

Figure  4-12  shows  the  monostatic  DD  term  comparison  at  3  GHz, 
where  the  effect  of  the  DD  term  is  more  apparent.  By  excluding  the  DD 
term,  the  divergence  of  the  solution  near  180°  is  alleviated,  but  the 
degradation  of  the  pattern  elsewhere  has  increased  greatly.  The  fact  that 
the  pattern  changes  when  the  DD  term  is  excluded  signifies  the 
importance  of  the  double  diffraction  mechanism  in  the  RCS  of  the  dihedral 
corner  reflector.  The  regions  where  the  DD  mechanism  is  more 
prominent  are  shown  by  larger  differences  between  the  DD  term  included 
and  excluded  patterns. 

Figure  4-13  shows  another  monostatic  DD  term  comparison  at  10 
GHz.  From  3  GHz  to  10  GHz  the  wavelength  shortens;  thus,  the  distance 
for  a  360°  phase  shift  also  shortens,  which  causes  the  RCS  pattern  to 
exhibit  more  oscillations.  The  increased  oscillation  makes  the  task  of 
pattern  comparisons  a  little  more  difficult.  However,  the  degradation  of 
the  RCS  pattern  from  excluding  the  DD  term  is  still  obvious,  especially  for 
theta  polarization.  All  of  the  monostatic  patterns,  regardless  of  frequency, 
have  the  same  problem  of  the  reflection  shadow  boundary  at  135°  and  180° 
causing  the  solution  to  diverge. 

What  happens  when  the  geometry  is  changed  to  a  bistatic  case?  The 
observation  point  leads  the  source  point  by  a  specified  angle  when  the 
bistatic  case  is  stated.  Figure  4-14  shows  the  90°  bistatic  RCS  of  the 
dihedral  at  3  GHz  for  theta  in  theta  out  and  phi  in  phi  out  polarization. 
With  the  DD  term  included  the  theta  in  theta  out  pattern  is  well  behaved, 
except  for  finite  discontinuities  at  180°  and  225°.  Excluding  the  DD  term 
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FIGURE  4-13.  MONOSTATIC  RCS  OF  9”  DIHEDRAL  AT  10  GHZ 
UTD  DOUBLE  DIFFRACTION  TERM  COMPARISON 


removes  the  abrupt  shifts  in  the  pattern,  but  degrades  the  remainder  of  the 
pattern.  With  the  DD  term  included  the  phi  in  phi  out  also  exhibits  finite 
discontinuities  at  180°  and  225°,  but  an  infinite  discontinuity  exists  at  270°. 
Here  the  source  is  at  180°  which  is  looking  down  the  center  of  the  dihedral 
and  the  observation  point  is  directly  to  the  side.  The  problem  arises  from 
the  fact  that  edge  1  is  on  the  reflection  shadow  boundary  of  edge  2,  and  the 
observation  point  is  also  on  the  RB  as  well  as  on  the  shadow  boundary  of 
edge  1,  as  seen  in  Figure  4-15.  Excluding  the  DD  term  does  not  alleviate 
the  problem,  and  even  worsens  the  divergence  of  the  solution  in  the  region 
near  270°. 


I 

x 

Figure  4-15  Bistatic  Reflection  Boundary 


Figure  4-16  shows  the  90°  bistatic  DD  term  comparison  at  10  GHz  for 
theta  in  theta  out  and  phi  in  phi  out  polarization.  Again,  the  finite 
discontinuity  is  present  at  225°  for  both  polarizations,  and  the  phi 
polarization  has  the  infinite  discontinuity  at  270°.  The  importance  of  the 
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double  diffraction  mechanism  for  the  bistatic  case  can  be  seen  by  the 
differences  between  the  RCS  patterns  when  the  DD  term  is  included  and 
excluded  at  both  3  GHz  and  10  GHz. 

No  improvement  is  gained  by  excluding  the  DD  term  for  the  bistatic 
case,  so  can  the  discontinuities  be  taken  care  of  another  way?  Possibly  the 
observer  could  be  taken  off  the  shadow  and  reflection  boundaries  by 
changing  the  bistatic  angle  to  an  angle  slightly  off  90°.  Figure  4-17  shows  a 
comparison  between  a  90°  and  89.7°  bistatic  angle  for  theta  in  theta  out  and 
phi  in  phi  out  polarization.  For  theta  polarization,  the  glitch  at  180°  was 
completely  removed,  and  the  glitch  at  225°  was  greatly  reduced,  all  without 
altering  the  overall  pattern.  For  phi  polarization,  the  glitch  at  180°  was 
worsened,  the  glitch  at  225°  was  reduced,  and  the  discontinuity  at  270°  was 
worsened  and  made  to  appear  like  the  90°  bistatic  without  the  DD  term. 

Shifting  the  bistatic  angle  a  little  more,  Figure  4-18  shows  a 
comparison  between  90°  and  89.5°  for  3  GHz.  For  theta  polarization,  the 
glitch  at  180°  was  removed,  the  glitch  at  225°  was  reduced,  but  now  a  new 
glitch  was  introduced  at  270°.  For  phi  polarization,  the  glitch  was  removed 
at  180°,  reduced  at  225°,  and  the  discontinuity  was  again  worsened  at  270°. 
Patterns  were  also  made  for  bistatic  angles  of  89.3°  and  89.0°,  but  they 
proved  to  be  too  far  off  of  90°  to  give  accurate  results,  and  are  not  included 
in  this  thesis.  The  89.7°  bistatic  yielded  best  results  for  theta  polarization, 
and  the  90°  bistatic  yielded  the  best  results  for  phi  polarization  at  3  GHz. 

Figure  4-19  shows  a  comparison  between  90°  and  89.7°  at  10  GHz. 
For  theta  polarization,  the  pattern  changed  most  at  225°  and  270°  when  the 
bistatic  angle  was  shifted  to  89.7°.  At  225°  it  is  difficult  to  tell  if  there  was 
pattern  improvement,  at  90°  the  glitch  was  in  the  positive  direction  and 
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FIGURE  4-18. 


BISTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 
90.0  AND  89.5  BISTATIC  ANGLE  COMPARISON 
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FIGURE  4-19 .  BISTATIC  RCS  OF  9"  DIHEDRAL  AT  10  GHZ 

90.0  AND  89.7  BISTATIC  ANGLE  COMPARISON 


when  the  bistatic  angle  changed  to  89.7°  the  glitch  reversed  to  the  negative 
direction.  Since  the  pattern  was  going  negative  before  225°  and  going 
positive  after  225°  leads  one  to  feel  more  confident  in  the  89.7°  pattern  since 
the  glitch  also  exhibits  the  same  characteristics  before  and  after  225°.  For 
phi  polarization,  the  10  GHz  patterns  responded  the  same  as  the  3  GHz 
patterns  when  the  bistatic  angle  was  changed  to  89.7°.  The  pattern  blew  up 
at  180°  and  the  glitch  was  removed  at  225°  when  the  bistatic  angle  was 
shifted  to  89.7°.  At  270°  the  discontinuity  was  not  reduced  by  the  new 
bistatic  angle,  instead  the  discontinuity  was  widened  and  appeared  like  the 
pattern  with  the  DD  term  excluded. 

Changing  the  bistatic  angle  to  89.5°,  Figure  4-20  shows  the 
comparison  with  90°  for  both  theta  and  phi  polarization.  For  theta 
polarization,  the  pattern  exhibited  the  same  change  at  225°  as  the  89.7° 
comparison,  but  elsewhere  the  pattern  differed.  The  discontinuity  at  270° 
was  unchanged,  and  minor  changes  appeared  throughout  the  pattern. 
The  extra  pattern  changes  that  were  not  experienced  by  the  89.7°  pattern 
and  the  patterns  at  3  GHz  were  again  due  to  the  fact  that  the  distance  for  a 
360°  phase  shift  is  shorter  at  10  GHz,  thus  more  sensitive  to  changes  in  the 
angles,  both  aspect  and  bistatic  .  For  phi  polarization,  the  pattern  was  well 
behaved  at  180°,  the  glitch  was  removed  at  225°,  and  the  discontinuity  was 
widened  at  270°.  The  phi  pattern  also  experienced  the  extra  pattern 
changes  due  to  the  angle  sensitivity.  The  best  patterns  for  the  bistatic  case 
at  10  GHz  were  the  89.7°  for  theta  polarization,  and  90°  for  phi  polarization. 
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FIGURE  4-20.  BISTATIC  RCS  OF  9"  DIHEDRAL  AT  10  GHZ 

90.0  AND  89.5  BISTATIC  ANGLE  COMPARISON 
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Now  that  the  MOM  and  UTD  codes  have  been  compared  with 
themselves,  how  do  they  compare  with  each  other?  The  MOM  is  more 
accurate  than  UTD,  but  the  increased  accuracy  of  MOM  has  the 
disadvantage  of  increased  computer  time.  The  UTD  code  is  two 
dimensional  so  it  cannot  provide  cross  polarization  information,  while  the 
MOM  provides  all  polarization  components.  At  2  GHz,  MOM  with  a  0.10 
wavelength  full  surface  patch  test  took  approximately  4  hours  of  CPU  time 
on  the  VAX  11/785  with  14  megabytes  of  main  memory,  whereas  the  UTD 
only  took  6  seconds  of  CPU  time.  One  run  with  the  MOM  code  provided 
patterns  for  all  polarization  combinations,  where  the  UTD  code  only 
provided  one  copolarization  combination. 

Figure  4-21  shows  the  comparison  between  MOM  and  UTD  for  the 
monostatic  RCS  at  2  GHz  for  theta  and  phi  polarization.  The  phi  patterns 
compared  very  well,  but  what  happened  to  the  theta  polarization?  The 
problem  was  that  the  UTD  code  only  calculated  the  diffraction  from  the 
outer  edges  of  the  dihedral,  and  did  not  include  diffraction  terms  from  the 
top  and  bottom  edges  since  there  are  no  top  and  bottom  edges  in  a  two 
dimensional  model.  With  the  MOM  UTD  comparison,  the  task  of  finding 
the  regions  where  these  missing  terms  were  most  needed  was  simplified. 
For  theta  polarization,  the  region  where  the  extra  terms  were  needed 
included  almost  the  entire  pattern.  This  was  reasonable  since  the  energy 
diffracted  when  the  electric  field  is  parallel  to  the  edge  is  less  than  when 
the  electric  field  is  perpendicular,  corresponding  to  theta  and  phi 
polarization.  This  was  exemplified  by  Figure  4-21,  where  the  UTD  theta 
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polarization  pattern  is  not  as  accurate,  and  the  UTD  phi  polarization  I 

pattern  is  accurate  when  compared  to  the  MOM  solutions.  I 

Figure  4-22  shows  the  comparison  between  MOM  and  UTD  for  the  1 

monostatic  case  at  3  GHz.  Again,  the  UTD  pattern  compares  well  with  the  ■ 

MOM  solution  for  phi  polarization  and  not  so  well  for  theta  polarization.  . 

Figure  4-23  shows  another  comparison  for  the  monostatic  case  at  10  GHz. 

The  UTD  departed  from  the  MOM  solution  more  for  the  theta  polarization 
than  the  phi  polarization,  but  the  phi  polarization  was  not  as  accurate  as  it 
was  for  2  GHz  and  3  GHz.  Normally  the  UTD  solution  should  improve  in 
accuracy  as  the  frequency  is  increased  since  UTD  is  a  high  frequency 
approximation.  However,  the  missing  diffractions  terms  had  a  greater 

effect  on  the  pattern  at  10  GHz  than  2  GHz  and  3  GHz.  \ 

Figure  4-24  shows  the  UTD  MOM  comparison  for  the  bistatic  case  at 
3  GHz.  The  UTD  did  not  compare  well  to  MOM  for  theta  polarization,  and 
did  not  compare  as  well  for  phi  polarization  as  the  monostatic  case.  The 
poor  comparison  indicated  that  more  diffraction  mechanisms  were  needed 
to  accurately  predict  the  RCS  at  bistatic  angles. 

Figure  4-25  shows  the  UTD  MOM  comparison  for  the  bistatic  case  at 
10  GHz.  Once  again  there  was  poor  correlation  between  the  patterns  for 
theta  polarization,  and  fair  correlation  for  phi  polarization.  The  two 
dimensional  UTD  proved  to  be  too  much  of  a  simplification  for  a  three 
dimensional  object  like  the  dihedral  corner  reflector. 
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FIGURE  4-23.  MONOSTATIC  RCS  OF  9"  DIHEDRAL  AT  10  GHZ 
MOM  VERSUS  UTD  COMPARISON 
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FIGURE  4-24.  90  DEGREE  BISTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 

MOM  VERSUS  UTD  COMPARISON 


Figure  4-26  shows  the  comparison  between  the  MOM  with  0.10 
wavelength  full  surface  test  patch  and  measurement  for  the  monostatic 
RCS  at  2  GHz  for  theta  and  phi  polarization.  For  theta  polarization,  the 
agreement  was  excellent  except  for  the  depth  of  the  null  at  110°,  which  was 
off  by  only  5  dBsm.  For  phi  polarization,  the  agreement  was  good,  but  not 
as  good  as  the  theta  polarization.  There  were  problems  with  the  null 
depths  at  22°,  67°,  and  82°,  but  the  difference  was  at  most  5  dBsm.  There 
were  other  areas  of  slight  discrepancy,  but  overall  the  pattern  compared 
well  for  phi  polarization.  Figure  4-27  shows  the  MOM  measurement 
comparison  for  the  theta  in  phi  out  cross  polarization  component.  The 
comparison  was  very  poor  between  measurement  and  the  MOM  for  the 
cross  polarization  component,  especially  when  compared  to  the  the 
copolarization  comparison.  One  explanation  for  the  poor  comparison  is 
the  measurement  system  was  approaching  the  limit  of  the  noise  floor  of 
the  receiver.  This  noise  is  shown  by  the  jaggedness  or  roughness  of  the 
measured  pattern.  Error  can  be  introduced  into  the  measurements  by  the 
cross  polarization  characteristics  of  the  transmitting  and  receiving 
antennas.  Coupling  between  the  transmitting  and  receiving  antennas  is 
also  a  problem  when  taking  monostatic  measurements.  Another  possible 
trouble  area  was  that  the  MOM  code  had  to  be  run  with  theta  equal  to  89°  to 
obtain  any  cross  polarization  information.  This  put  the  source-receiver  1° 
(smallest  angular  increment  allowed  by  the  ESP  sode)  off  the  principal 
plane.  When  the  code  was  run  at  theta  equal  90°,  the  cross  polarization 
component  was  absent,  so  the  code  was  rerun  with  the  small  change  in 
theta  and  then  the  cross  polarization  appeared.  Assuming  that  the 
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FIGURE  27.  MONOSTATIC  RCS  OF  9”  DIHEDRAL  AT  2  GHZ 
MOM  VERSUS  MEASURED  COMPARISON 


measurement  system's  target  mount  had  a  small  error  in  the  position  of 
the  dihedral,  running  the  code  with  the  89°  angle  was  not  a  bad 
approximation.  The  relative  level  of  the  MOM  and  measured  was  not  bad, 
just  the  actual  pattern  oscillations. 

Figure  4-28  shows  the  MOM,  with  0.25  wavelength  full  surface  test 
patch,  meaf  rement  comparison  for  the  monostatic  RCS  at  1 0  GHz  for 
theta  polarization.  The  comparison  is  excellent!  The  phi  polarization 
pattern  was  not  measured  by  OSU-ESL,  but  the  MOM  solution  is  shown  for 
completeness.  Figure  4-29  shews  the  cross  polarization  comparison  for 
theta  in  phi  out  polarization.  Both  measurements  were  for  theta  in  phi  out 
polarization,  and  the  fact  that  the  patterns  were  not  identical  points  to 
problems  with  the  measurement  process.  The  MOM  code  was  run  again 
with  theta  equal  to  89°  to  obtain  cross  polarization  information.  Figure  4-30 
shows  the  phi  in  theta  out  polarization  MOM  solution,  for  which  there  was 
no  measurement  taken  for  comparison.  For  the  monostatic  case,  the  theta 
in  phi  out  polarization  pattern  is  identical  to  phi  in  theta  out  polarization 
pattern. 

Figure  4-31  shows  the  MOM,  with  0.15  wavelength  full  surface  patch 
test,  measurement  comparison  for  the  90°  bistaU.  case  at  3  GHz  for  theta 
and  phi  polarization.  The  comparison  for  theta  polarization  is  excellent! 
The  phi  polarization  pattern  agreement  is  good,  but  not  as  good  as  the  theta 
comparison  in  terms  of  predicting  the  depth  of  the  nulls.  The  phi  pattern 
again  has  problems  in  the  depth  of  the  nulls  at  142°,  162°,  200°,  210°,  and 
282°.  Figure  4-32  shows  the  cross  polarization  comparison  for  theta  in  phi 
out  and  phi  in  theta  out  polarization.  The  comparison  is  not  as  good  as  the 
copolarization  comparisons  were.  For  theta  polarization,  the  patterns 
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compare  well  until  about  90°  where  they  deviate  until  21 0°  when  they  come 
together  again.  For  phi  polarization,  the  pattern  shapes  compare  okay,  but 
for  some  reason  the  magnitude  of  the  MOM  solution  is  5  to  15  dBsm  higher 
throughout  the  pattern.  The  patterns  don’t  deviate  much  until  the  region 
beyond  330°. 

Figure  4-33  shows  the  MOM  measurement  comparison  for  the  90° 
bistatic  case  at  10  GHz  for  theta  and  phi  polarization.  The  MOM  theta 
polarization  pattern  compared  well  with  the  measured  pattern,  except  for 
minor  differences  mostly  in  the  region  between  215°  and  245°.  The  MOM 
and  measured  patterns  for  phi  polarization  also  compared  well,  except  for 
the  region  near  290°  to  295°.  Figure  4-34  shows  the  cross  polarization 
comparison  between  MOM  and  measured  at  10  GHz.  For  theta  in  phi  out 
polarization  the  peaks  were  accurately  predicted,  but  the  remainder  of  the 
pattern  didn't  compare  as  well.  For  phi  in  theta  out  polarization  the  peaks 
weren't  predicted  as  accurately  as  for  theta  in  phi  out  polarization,  but  the 
rest  of  the  pattern  compared  better.  The  cross  polarization  comparisons 
for  the  bistatic  case  at  10  GHz  were  the  best  out  of  all  of  the  cross 
polarization  comparisons  for  MOM  versus  measured. 
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The  Electromagnetic  Surface  Patch  code  predicted  all  polarization 
components  accurately  for  the  dihedral  comer  reflector  for  the  various 
frequencies  and  geometries  of  the  measurements.  Confidence  in  the  ESP 
code  allows  other  geometries  for  the  dihedral  to  be  explored.  Up  to  this 
point  all  models  had  been  run  in  the  theta  equal  89°  or  90°  elevation  plane. 
What  happens  to  the  RCS  if  the  source  and  receiver  is  moved  off  the 
principal  plane?  Figure  4-35  shows  the  RCS  for  the  monostatic  case  at  3 
GHz  for  theta  and  phi  polarization  for  various  elevation  angles.  Changing 
theta  from  89°  to  75°  greatly  reduced  the  backscatter  for  both  theta  and  phi 
polarization.  However,  moving  further  off  the  principal  plane  did  not 
reduce  the  copolarization  backscatter  very  much.  Figure  4-36  shows  the 
cross  polarization  component  of  the  RCS  of  the  dihedral  at  3  GHz  for  the 
monostatic  geometry.  The  cross  polarization  component  was  lowest  for  the 
theta  equal  89°  pattern.  When  moving  of  the  principal  plane  to  theta  equal 
75°,  60°,  and  45°,  the  cross  polarization  component  increased.  Also  shown 
in  Figure  4-36  is  an  RCS  pattern  for  the  azimuth  angle  fixed  at  145°  and  the 
variable  angle  is  elevation  (theta).  The  pattern  comparison  shows  both 
copolarization  and  cross  polarization  for  the  dihedral  at  3  GHz  for  the 
monostatic  case.  For  the  most  part  the  copolarization  components  were 
higher  in  magnitude  than  the  cross  polarization  components,  but  in  some 
regions  the  cross  polarization  component  was  as  large  and  sometimes 
larger  than  the  copolarization  component. 
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FIGURE  4-35.  MONOSTATIC  RCS  OF  9"  DIHEDRAL  AT  3  GHZ 
OFF  PRINCIPAL  PLANE  PATTERN  COMPARISON 


The  UTD  model  employed  in  this  thesis  was  more  efficient  in  terms 
of  computer  time  and  resource  consumption.  However,  the  UTD  model 
was  less  accurate  due  to  its  two  dimensional  limitation.  The  two 
dimensional  UTD  model  didn’t  accurately  predict  the  RCS  of  the  dihedral 
corner  reflector  because  the  dihedral  is  a  three  dimensional  object  with 
finite  edges  and  corners.  The  UTD  model  ignored  some  of  the  diffraction 
effects  which  are  present  in  a  three  dimensional  object. 

The  three  dimensional  MOM  model  implemented  was  more 
accurate  for  predicting  the  RCS  of  the  dihedral  comer  reflector.  But  the 
accuracy  of  the  MOM  model  was  offset  by  an  increased  cost  of  computer 
time  and  resource  utilization.  The  MOM  provided  cross  polarization  RCS 
information,  whereas  the  UTD  model  couldn’t  predict  the  cross 
polarization  RCS  because  of  its  two  dimensional  nature.  The  MOM  model 
had  an  upper  limit  in  the  electrical  size  of  the  object  it  could  analyze, 
translating  into  a  maximum  frequency  and  physical  size  of  the  object. 
Above  the  maximum  frequency  and  physical  size  of  the  object,  the  MOM 
model  required  more  resources  than  the  computer  had  to  offer. 

The  copolarization  RCS  measurements  taken  at  The  Ohio  State 
University  ElectroScience  Laboratory  (OSU-ESL)  compared  well  with  the 
MOM  model's  results.  However,  the  cross  polarization  RCS 
measurements  didn’t  match  up  with  the  MOM  model  as  well  as  the 
copolarization  RCS  measurements.  The  cross  polarization  RCS 


measurements  were  not  as  smooth  and  continuous  as  expected;  instead 
they  exhibited  an  irregularity  which  brought  doubts  as  to  their  accuracy. 
The  magnitudes  of  the  cross  polarization  RCS  measurements  weren't  too 
far  above  the  stated  noise  floor  of  the  OSU-ESL  compact  range.  A  problem 
with  cross  polarization  RCS  measurements  is  the  depolarization 
characteristics  of  the  transmitting  and  receiving  antennas.  Coupling  is 
another  problem  with  the  antennas  when  they're  in  close  proximity  to  each 
other,  as  in  the  monostatic  measurements.  The  target  mount  which  holds 
the  dihedral  corner  reflector  can  also  introduce  error  in  the  measurement 
due  to  its  depolarization  characteristics,  especially  for  bistatic  angles. 

The  MOM  model  predicted  a  substantial  cross  polarization  RCS  for 
the  dihedral  corner  reflector  when  the  source  and  observation  point  moved 
away  from  the  principal  plane.  The  dihedral  corner  reflector  would 
normally  be  viewed  away  from  the  principal  plane  when  used  as  a 
component  on  a  complex  target  model.  No  off  principal  plane 
measurements  were  taken  for  comparison,  but  the  accuracy  of  the  MOM 
model  has  been  verified.  A  significant  cross  polarization  RCS  is  desired  so 
that  this  information  can  be  used  as  an  additional  parameter  in  target 
identification  in  future  radar  systems. 

Recommendations 

First,  the  two  dimensional  UTD  model  should  be  extended  to  the 
three  dimensional  problem  for  the  reason  that  the  UTD  model  is  much 
more  efficient  in  terms  of  computation  time.  This  extension  would  be  in 
the  area  of  additional  diffraction  and  reflection  mechanisms  not 
incorporated  in  the  two  dimensional  model.  Edge  diffraction  mechanisms 
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need  to  be  included  for  the  top  and  bottom  edges  of  the  dihedral  comer 
reflector,  whereas  the  two  dimensional  model  only  included  diffraction 
mechanisms  for  the  two  outer  edges.  Comer  diffraction  mechanisms  need 
to  be  included  if  accurate  cross  polarization  information  is  desired. 

Second,  the  MOM  model  should  be  modified  to  increase  the 
computational  efficiency.  Different  basis  (expansion)  and  testing  (weight) 
functions  might  reduce  the  computer  time  required  to  compute  the 
impedance  and  excitation  matrices.  Entire  domain  basis  functions  could 
possibly  reduce  the  computation  time.  Symmetry  within  the  impedance 
matrix  might  exist,  which  could  be  exploited  to  reduce  the  number  of 
elements  to  compute.  Another  alternative  would  include  the  use  of  a  faster 
computer  that  has  a  large  core  (main)  memory.  However,  modifying  the 
MOM  model  should  be  investigated  before  brute  force  is  implemented. 

Third,  other  computational  techniques  for  electromagnetic 
scattering  should  be  explored.  The  Spectral  Domain  and  Conjugate 
Gradient  techniques  have  been  identified  as  possible  alternatives  to  the 
UTD  and  MOM  techniques.  The  first  area  to  examine  is  the  application  of 
the  technique  to  open  body  surfaces  like  the  dihedral  comer  reflector. 
Many  techniques  are  limited  to  closed  body  surfaces,  such  as  cylinders  and 
cones. 

Finally,  other  structures  should  be  modeled  for  cross  polarization 
RCS  for  both  monostatic  and  bistatic  geometries  off  the  principal  plane.  A 
flat  plate,  trihedral  corner  reflector,  cone,  cylinder,  and  ogive  are 
examples  of  other  structures  to  model.  If  the  object  is  a  closed  body  surface 
then  other  computational  techniques  become  available  to  the  researcher, 
such  as  the  MOM  solution  of  the  Magnetic  Field  Integral  Equation  (MFIE). 


86 


Summary 


This  thesis  examined  the  Radar  Cross  Section  (RCS)  of  the  dihedral 
corner  reflector  for  both  monostatic  and  bistatic  geometries.  The  full 
polarization  of  the  RCS  was  maintained  while  the  dihedral  corner  reflector 
was  modeled  on  and  off  the  principal  (horizontal  reference)  plane.  The 
dihedral  corner  reflector  was  composed  of  two  perfectly  conducting  9  inch 
square  flat  plates  joined  along  a  common  edge  so  that  they  formed  an 
interior  angle  of  90  degrees.  The  frequencies  of  interest  included  2,  3,  and 
10  GHz.  The  Uniform  Theory  of  Diffraction  and  the  Method  of  Moments 
were  utilized  in  the  modeling.  The  results  from  the  models  were 
compared  with  measurements  for  the  principal  plane  RCS.  The  MOM 
model  was  then  extended  to  off  principal  plane  RCS,  which  was  not 
measured. 
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Appendix  A:  Definition  of  Line  Types  and  Terms 
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Figure  A-l.  Definition  of  Line  Types 
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Definition  of  Terms 


£ 

H 

1 

s. 

o 

A 

h 

A 

V 

I 

A 

D 

<He(6J)> 

He 

<x,y> 


Complex  Electric  Field  Intensity  Vector 
Complex  Magnetic  Field  Intensity  Vector 

Complex  Electric  Current  Density  Vector 
Complex  Scattering  Coefficient 

Radar  Cross  Section 
Horizontal  Unit  Vector 
Vertical  Unit  Vector 

Incident  Propagation  Direction  Unit  Vector 

Diffraction  Propagation  Direction  Unit  Vector 
Diffraction  Dyadic 

Diffraction  Coefficient 

Inner  Product 
Weight  or  Test  Function 

Expansion  or  Basis  Function 
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Appendix  B:  Measurement  Facility 


The  measurements  were  performed  at  the  The  Ohio  State  University 
ElectroScience  Laboratory  (OSU-ESL)  compact  scattering  range.  The 
compact  range  is  a  broadband,  coherent,  stepped-CW  radar  measurement 
facility  which  has  the  capability  of  measuring  either  \,he  complex 
scattering  spectrum  of  a  target  over  a  broad  range  of  frequencies  at  a  fixed 
look  angle,  or  of  measuring  the  angular  RCS  pattern  at  the  fixed 
frequency.  The  compact  range  has  the  ability  to  measure  both  amplitude 
and  absolute  phase  over  an  extremely  broad  band  of  1  to  18  GHz. 

In  either  the  frequency  scan  or  angular  scan  mode,  each  data 
sample  is  distilled  from  a  set  of  four  separate  measurements: 

1)  target 

2)  target  background 

3)  calibration  target 

4)  chamber  background 

For  a  target  in  free  space,  2  and  4  are  the  same;  for  a  target  mounted  on  a 
ground  plane,  2  is  a  measurement  of  the  scattering  of  the  ground  plane 
alone.  An  array  of  these  four  measurements  is  stored  in  the  computer  and 
a  post  facto  vector  subtraction  is  used  to  remove  undesired  deterministic 
signals  due  to  the  chamber,  ground  plane,  horn  coupling,  pedestal,  etc., 
from  the  desired  target  signal.  At  each  sample,  the  result  is  calibrated 
against  the  known  RCS  of  the  calibration  object,  usually  a  sphere.  The 
procedure  has  been  shown  to  produce  excellent  results  for  free  space 
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measurements  of  a  variety  of  canonical  objects.  When  the  measured 
results  are  compared  to  the  exact  solutions  for  these  objects,  the  error 
component  is  some  30  dB  below  the  target  response. 

The  measurement  system  is  shown  in  block  form  in  Figure  B-l. 
Data  acquisition  is  controlled  by  a  VAX  11/750  minicomputer.  The  system 
presently  is  able  to  acquire  approximately  4  samples  per  second  in  the 
frequency  scan  mode  and  8  samples  per  second  in  the  angular  scan  mode. 

The  source  is  a  Watkins-Johnson  crystal  stabilized  frequency 
synthesizer  with  a  transmit  power  capability  of  approximately  100 
milliwatts  (mW).  A  second  Watkins-Johnson  is  used  to  downconvert  the 
received  signal  to  1  GHz,  bypassing  the  down-conversion  loop  in  the 
Scientific  Atlanta  receiver  and  stabilizing  the  IF  noise  level.  The  resulting 
receiver  sensitivity  is  approximately  -100  dBm  over  the  entire  band  from  2 
to  18  GHz. 

The  antenna  feed  system  uses  a  pair  of  broadband,  coaxial,  TEM 
horns  in  a  side-by-side  configuration  at  the  focal  point  of  a  12’  diameter 
Scientific  Atlanta  offset  reflector  to  which  4’  flared  edges  designed  at  OSU- 
ESL  have  been  attached  in  order  to  eliminate  the  interference  of  edge 
effects  with  the  plane  wave  field.  The  plane  wave  illuminated  sweet  spot  is 
8'  wide  by  4’  high. 

The  target  is  supported  on  a  low  cross  section  metallic  pedestal 
designed  at  OSU-ESL  with  a  500  pound  capacity  and  computer  controlled 
rotation  capability  accurate  to  within  0.01  degrees.  The  anechoic  chamber 
is  40  by  20  by  60  feet  with  ridged  absorber  on  the  floor  and  ceiling,  and 
pyramidal  cone  panels  on  the  walls. 
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In  a  typical  frequency  scan,  the  sampling  rate  is  every  10  MHz  so 
that  the  total  open  time  window  is  100  nanoseconds.  With  the  target 
located  some  30'  from  the  back  wall  and  24’  from  the  main  dish,  this 
sampling  rate  avoids  aliasing  of  the  chamber  return  inti*  th«  target  area. 
In  the  most  recent  version  of  the  radar  system,  the  stepped-CW  mode  has 
been  replaced  by  a  pulsed  stepped-CW  mode  of  OSU-ESL  design  which  uses 
a  set  of  switches  to  open  and  close  the  transmitter  and  receiver  with  a  40 
nanosecond  wide  pulse.  With  this  system,  the  return  from  the  back  wall, 
main  reflector,  and  direct  hom-to-hom  coupling  is  eliminated  entirely  so 
that  the  system  can  devote  nearly  all  of  the  -80  dB  dynamic  range  available 
in  the  receiver  to  the  target  response.  In  addition,  the  pulsed  system 
allows  one  rather  than  two  antennas  to  be  used  at  the  feed,  thus  improving 
the  focus  of  the  offset  fed  reflector  configuration. 

The  AEL  feed  antennas  operate  from  1  to  12  GHz  or  2  to  18  GHz  in  a 
single  frequency  scan.  A  typical  spectral  bandwidth  of  2  to  18  GHz  gives  a 
time  resolution  of  0.0625  nanoseconds  on  a  spatial  resolution  of  0.375 
inches  downrange.  A  resolution  of  this  magnitude  allows  the  finer  details 
of  the  target’s  structural  scattering  characteristics  to  be  isolated  into 
distinct  local  scattering  centers  and  diagnosed  separately. 

The  sensitivity  of  the  measurements  is  ultimately  limited  not  by  the 
noise  floor  of  the  receiver  but  by  the  clutter  floor,  which  is  determined  by 
the  ability  to  achieve  the  vector  subtraction.  Since  the  chamber  clutter 
from  points  at  large  distances  from  the  target  region  may  be  removed 
either  numerically  or  electronically  with  the  pulsed  system,  it  is  the  clutter 
which  returns  in  or  near  the  target  time  frame  which  sets  the  clutter  floor. 
As  stated  above,  for  free  space  measurements  this  error  has  been  found  to 
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be  approximately  30  dB  below  the  target  response.  When  the  target  is 
mounted  on  a  ground  plane,  it  is  the  ground  plane  scatter  which  primarily 
determines  the  clutter  floor.  The  subtractability,  not  just  the  absolute  level, 
of  the  ground  plane  return  is  the  key.  Therefore,  the  stability  and 
positional  repeatability  of  the  ground  plane  are  the  limiting  factors.  In  the 


past  year,  ground  plane  clutter  levels  of  -55  dBsm  absolute  have  been 


achieved  routinely. 
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FIGURE  B-l.  COMPACT  RANGE  BLOCK  DIAGRAM 
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